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Abstract 


method  of  coupled  azimuthal  potentials  (CAP)  was  applied  to  a 
waveguide  model  of  an  zueially  symmetric  engine  inlet  to  cuialyze  the  fields 
in  the  region  where  the  front  face  of  the  engine  terminates  the  waveguide. 
Appropriate  boundary  conditions  were  derived  and  the  finite  element  method 
was  used  to  solve  for  the  potentials. 

Lagr^ulgian  of  the  CAP  equations  does  not  provide  for  the  enforce¬ 
ment  of  Neunumn  boundary  conditions.  This  prevents  exact  implementation 
of  the  correct  boundzury  conditions  for  the  azimuthal  magnetic  field. 
Dirichlet  boundary  conditions  for  the  azimuthal  electric  and  maignetlc 
fields  were  enforced  for  a  standing  wave  condition  in  the  inlet  model 
with  a  conducting  flat  plate  termination.  The  computed  values  for  the 
Interior  field  components  were  compared  by  evaluating  the  standard 
deviation.  Three  trials  were  performed  with  varying  finite  element  mesh 
densities.  It  was  found  that  as  the  mesh  density  Increased,  the  stand¬ 
ard  deviation  for  the  computed  field  components  decreased. 


A 

from  an  interpolation  of  the  error  measurements,  it  was  determined  that 

an  extremely  fine  mesh  is  required  for  eui  acceptable  error  (approx- 

2 

Imately  7500  trlangles/X^  for  2%  error) .  For  this  reason  the  finite 
element  implementation  of  the  CAP  equations  for  this  jet  engine  inlet 
model  is  unfeasible. 


A  FINITE  ELEMENT  COMPUTATION  OF  THE 


ELECTROMAGNETIC  FIELDS  WITHIN  AN 
ENGINE  INLET  MODEL 

I .  Introduction 


Background 

The  jet  engine  Inlet  provides  a  significant  contribution  to  the 
overall  radar  cross  section  of  a  modern  aircraft.  Before  this  contribu¬ 
tion  can  be  effectively  reduced,  the  electromagnetic  wave  scattering 
phenomena  from  this  region  must  be  accurately  evaluated  euid  understood. 
This  is  a  relatively  difficult  problem  due  to  the  generally  complex 
nature  of  the  geometry  of  the  inlet.  However,  if  an  understzuiding  of 
the  scattering  mechanisms  from  a  simple  intake  geometry  is  developed, 
khe  insight  could  be  auialytlcally  applied  to  the  complex  geometry  of  an 
actual  inlet. 

The  engine  Intake  usually  consists  of  a  hollow  tubular  structure 
terminatea  by  the  fan  blades  of  the  eitgine  compressor.  It  is  well  known 
that  these  fan  blades  significantly  affect  the  radar  return  by  modulating 
the  echo  (1:435) .  For  this  reason,  any  simplification  in  the  geometry 
of  an  inlet  model,  would  have  to  provide  for  the  bladed  structure  of  the 
compressor.  The  model  to  be  used  here  will  consist  of  a  hollow  circular 
cylindrical  waveguide  terminated  by  a  conducting  cone  centered  on  a  flat 
plate.  Boundary  conditions  for  the  flat  plate  will  be  appropriately 
chosen  to  simulate  the  fan  blade  configuration. 

Problem 

The  basic  problem  is  to  evaluate  the  electromagnetic  wave  inter- 
aotimis  with  a  simplified  termination  model  of  a  jet  engine  inlet. 


Ihla  study  will  attempt  to  determine  the  scattered  electromagnetic 
fields  from  a  proposed  inteUce  termination  model.  No  attempt  will  be  made 
to  recommend  methods  of  reducing  the  radeu:  cross  section  of  a  jet  engine. 

Literature  Review 

Several  authors  have  previously  developed  and  analyzed  models  of  the 
jet  engine  Intedce  area  in  an  attempt  to  provide  insight  into  the  electro- 
sugnetlc  scattering. 

J<^  W.  Moll  emd  Rolf  6.  Seecamp  aneJ.yzed  an  inlet  modelled  by  a 
terminated  circular  cylindrical  waveguide.  The  termination,  which  repre¬ 
sents  the  Initial  stages  of  a  compressor,  consists  of  a  pair  of  conducting 
planar  blade  structures  as  depicted  in  Figure  1.  Each  structure  has  a 
different  number  of  blades  and  blade  widths.  The  field  incident  on  the 
opcmlng  of  the  inlet  is  represented  by  a  sum  of  inwardly  travelling  TE 
nodes,  ihe  backscattered  field  from  the  termination  is  expressed  as  a 
sum  of  TE  modes  travelling  towiurd  the  opening  of  the  duct.  Any  TH  modes 
excited  by  the  termination  are  neglected.  Along  the  surface  of  the 
blades,  the  tangential  electric  field  was  set  to  zero.  Using  this 
boundaury  condition,  a  set  of  equations  was  obtained  to  solve  for  the 
scattering  coefficients  for  the  modes  generated  at  the  terminatlm.  The 
total  field  at  the  opening  of  the  inlet  is  then  equated  to  the  sum  of 
the  incident  and  reflected  modes.  The  Stratton-Chn  integral  is  then  ap¬ 
plied  to  obtain  the  scattered  field  from  any  aspect  angle  (2) . 

A  q^lte  similar  model  was  analyzed  by  P.  H.  Pathak  and  C.  C.  Huang. 

It  consists  of  a  hollow  circular  cylindrical  waveguide  terminated  by  a 
single  pla;  r  blade  ructure  as  in  Figure  2.  All  of  the  blades  are 


Figure  2.  Inlet  Model  Used  by  Pathak  and  Huang 


assumed  to  be  of  the  same  width  separated  by  a  finite  cuigular  distance. 
Basically  their  approach  is  to  calculate  scattering  matrices  at  the 
opening  of  the  duct  (Sj.!*  ^22  ^21^  ^  both  an  externally 

generated  Incident  field  and  the  reflected  field  due  to  the  blade 
structure  discontinuity.  All  of  these  matrices  aze  calculated  using 
the  Uniform  Theory  of  Diffraction.  The  scattering  matrix  at  the 
termination  is  evaluated  by  a  geometrical  optics  approximation  to  the 
current  distribution  (Jg  =  2«d  x  H)  on  the  fan  blades.  The  radcur  cross 
section  of  this  model  is  then  determined  via  the  Multiple  Scattering 
Method,  which  relates  the  backscattered  field  to  the  Incident  field  (3) . 

T.W.  Johnson  and  D.  L.  Moffatt  analyzed  the  wave  scattering  from  a 
circular  cylindrical  waveguide  with  a  non-plan2u:  termination  as  illus¬ 
trated  in  Figure  3.  The  termination  considered  is  an  axially  symmetric 
cone  centered  on  a  flat  plate.  Both  the  cone  and  plate  are  perfectly 
conducting.  Using  the  Wiener-Hopf  technique,  the  scattering  matrices 
for  the  opening  of  the  waveguide  are  calculated.  The  dyadic  magnetic 
Green's  function  and  the  Method  of  Moments  are  then  used  to  numerically 
compute  the  current  distribution  on  the  cone.  From  this  distribution, 
the  reflection  coefficient  for  the  cone  was  determined  for  the  TBn  mode 
The  radar  cross  section  is  then  evaluated  usiiig  this  backscatter  coef¬ 
ficient  and  the  scattering  matrices  calculated  for  the  open  end  of  the 
waveguide  (4) . 

Approach 

The  inlet  model  that  will  be  analyzed  represents  a  slight  modifica¬ 
tion  of  the  model  used  by  Johnson  and  Moffat.  It  consists  of  a  circular 
cylindrical  waveguide  terminated  by  a  perfectly  conducting  cone  centered 
on  a  flat  plate.  However,  the  boundeury  conditions  enforced  on  the  plate 


Hon-planar  Temlxiatlon  Used  by  Johnson  and  Moffat 


will  not  be  those  of  a  perfect  conductor,  instead,  a  set  of  boundary 
conditions  will  be  used  to  make  the  flat  plate  resemble  the  blade 
structure  of  the  compressor  termination  in  cm  actual  inlet. 

Since  fan  blades  eure  separated  by  a  finite  anguleu:  distance,  emy 
Current  induced  on  the  blades  cannot  flow  azimuthally.  To  simulate  this 
effect  in  the  model,  the  radially  directed  tangential  magnetic  field  at 
the  surface  of  the  flat  plate  will  be  set  to  zero.  'Riis  will  eliminate 
azlfflutlud.  ctirrent  in  the  plate  and  make  it  appeeu:  as  a  set  of  infinitely 
thin  wires  directed  radially  outward.  The  radially  directed  tcmgential 
electric  field  at  the  surface  of  the  plate  will  also  be  set  to  zero, 
which  is  the  appropriate  boundary  condition  for  a  conductor. 

With  the  aid  of  the  Coupled  Azimuthal  Potential  (CAP)  formulation, 
the  finite  element  method  emd  the  unlmoment  method,  the  total  fields 
within  the  model  will  be  determined.  The  CAP  consists  of  two  second 
order  differential  equations  developed  by  Morgan,  Chang  and  Med  to  solve 
scattering  problems  involving  generally  lossy  isotropic  Inhomogenous 
rotationally  symmetric  media  (5) .  These  potentials,  which  eu:e  directly 
related  to  the  electric  and  magnetic  fields,  will  be  numerically  computed 
using  the  finite  element  method.  This  requires  that  an  Integral  func¬ 
tional  be  used  to  set  up  the  numerical  solution  for  the  fields  in  the 
waveguide.  The  functional,  when  minimized,  is  analytically  equivalent 
to  enforcing  the  differential  equations.  Finally,  the  unlmoment  method 
will  be  used  to  determine  a  linear  combination  of  fields  which  will 
satisfy  the  boundary  conditions  of  the  inlet  model. 

Chapter  II  will  be  devoted  to  presenting  the  CAP  formulation  theory 
and  the  finite  element  implementation  of  the  differential  equations.  In 
addition,  the  application  of  the  unimoment  method  to  this  problem  will 
be  detailed. 
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Chapter  III  will  give  the  results  ^md  Chapter  IV  the  conclusions 


2uid  recommendations  resulting  from  the  analysis  of  the  inlet  model. 


Analytical  Formulation 


This  chapter  will  discuss  the  analysis  used  in  obtaining  the  system 
of  equations  for  evaluation  of  the  electromagnetic  fields  in  the  interior 
of  the  engine  inlet  model.  The  derivation  of  the  coupled  azimuthal  po¬ 
tential  (CAP)  formulation  will  be  presented.  Numerical  computation  of 
the  CAP  by  use  of  the  finite  element  method  will  also  be  developed.  In 
addition,  the  application  of  the  unimoment  method  to  calculate  the  total 
fields  will  be  discussed. 


CAP  Formulation 

The  CAP  is  valid  in  generally  inhomogenous  isotropic  rotationally 
symmetric  media  (bodies  of  revolution)  as  shown  in  Figure  4.  Its  basic 
utility  is  that  it  converts  a  three-dimensional  scattering  problem  into 
one  involving  only  two  dimensions.  Since  the  medi\m>  is  axl-syimietric, 
the  fields  have  a  known  functional  behavior  with  respect  to  the  ^  coord¬ 
inate  (in  a  standard  cylindrical  coordinate  system) .  Therefore,  if  the 
fields  are  determined  for  a  constemt  R-Z  cross  section,  they  are  known 
everywhere  in  the  media. 


Z 


Figure  4.  Meridional  Two-dimensional  Solution  Domain  (6:203) 


The  CAP  formulation  development  begins  by  defining  a  Cartesian 
coordinate  system  consisting  of  normalized  circular  cylindrical  co¬ 
ordinates.  All  points  within  the  system  may  be  given  by  (R,Z,(^)  where 
R  «  XoPf  Z  s  KqZ,  ^  ^  and  Kq  »  2ir/Xo  is  the  free  space  wave  number 

of  the  fields.  The  fields  are  then  exp2mded  into  cut  exponential 
Fourier  series  in  the  ^  emrdlnate: 

E(K,Z,»  -  exp  (Jn*) 

n  ;(R,z,«  -  ? 

o  ** 

where  *  y^/e^.  When  these  expeutsions  are  placed  in  Maxwell's  equa¬ 

tions,  V  X  E  •  -j(oyH  and  V  x  H  »  jotcE,  the  following  coupled  system  of 
equations  is  obtained: 


(1) 

(2) 


^  [ 


3h2,n  _  3hR 
3R  3Z 


,  h  -  j  riSRJL-  isgiill 

^r"^,n  ^  I  iz  .3R  J 


r  3(Rhd.n)  _ 


»<'r«R,n-  3 


"■A.n  -3  [j-z.n  - 
[j"Vn  - 


“r'z,n  -  3 


RM, 


(3) 

(4) 

(5) 

(6) 

(7) 

(8) 


Substituting  among  those  equations,  it  can  easily  be  shown  that 


(9) 


e^(R,2)  »  jfjj(nVi|»2^  +  X  7l|»2)  +  ^(ipj^/R) 

h^(R,Z)  »  jfjj(nVt()2  -  Re^^  x  +  (f(i|i2/R)  <10) 

irtiere  ^>./R  *  e.  ,  ^./R  "  h.  and  the  gradient  operator,  V,  is  defined  as 

X  “  Yf^ 

A  A 

R  3/3R  +  Z  3/3Z.  Addltlomaiy, 

f^(R»Z)  -  u,(R,Z)  e^(R,Z)R^  -  T?  (11) 

From  Eqs  (9)  and  (10)  the  radial  and  Z  directed  components  of  the  elec¬ 

tric  emd  magnetic  field  may  be  determined.  These  components  may  be 
shown  to  be 


(12) 

.,(*,«  -  ^  -  R 

(13) 

V®'®'  *  ^*n‘®  ^  ^ 

(14) 

hj(R,«  -  Jt„(^  +  R  ^ 

(15) 

The  coupled  potentials,  and  satisfy  a  system  of  peurtlal  dif¬ 
ferential  equations  given  by: 

A 

i  fj^CRCr^^’l  +  X  7<i2^  ® 

A 

V  *u(RMjf^*2  “  “♦  *  ^*1^  ♦'r  *2^  “  ® 

To  utilize  the  finite  element  method  to  numerically  compute  the  potentials, 
it  is  necessary  to  develop  a  Lagrangian  functional.  The  stationary  point 
Of  the  Integral  of  this  functional  is  analytically  equivalent  to  the 
partial  differential  equations.  For  this  set  of  differential  equations 


II 


L  ■  f  V(|»,  (Re  7ij»  +  n«>  X  74)  )  + 
n  X  X  j.  i 


74)^  *  (Ry^  74)^  -  n4»  X  74)^^)  -  (14) 


Finite  Element  Method 

The  finite  element  method  is  a  numerical  method  of  calculating  the 
approximate  solution  of  a  partial  differential  equation  in  em  enclosed 
area.  It  does  so  by  assuming  a  form  of  the  solution  (linear,  quadratic, 
etc.)  within  a  set  of  "finite  elements"  which  span  the  total  eurea.  By 
integrating  am  Integral  functional  over  each  element,  the  displacement  of 
a  single  nodal  value  on  em  element  is  related  to  adjacent  nodal  values. 
This  process  results  in  a  system  of  equations.  The  simultaneous  solution 
of  this  set  of  equations  results  in  the  values  of  nodal  coefficients, 
thus  providing  a  plecew|.se  approximation  of  the  solution  from  node  to 
node. 

2  2 

When  utilizing  the  Lagranglan,  the  term  ***  produces 

a  singularity  €ilang  the  R  »  0  eucls.  To  overcome  this  difficulty,  a 

change  of  veurieUsle  must  be  performed  on  the  scalar  potentials,  letting 

Re.  «  4),  and  Rh^  *■  4’.>«  It  then  becomes  necessary  that  the  boundary 
f,n  1  9,n  2 

conditions  satisfied  by  the  azimuthal  field  quantities  along  R  •  0  be 
Specified.  These  bound^u7  conditions  are: 


R^  -  Ir-0  “  0  for  n  ^  1 


12£ijl  +  (  3i»Sy.)  e  1 
DR  '  3R  '  ♦,n  I 


0  for  n  -1 


/  iiSHr.)  h.  I 


■■  0  for  n  >1 


.2 


2 

The  boundeury  conditions  for  n  «  1  reduce  to  homogenous  Neum^uln  boundary 
conditions  when  the  relative  permittivity  and  permeability  are  invariant 
with  respect  to  R  (5:416).  To  enforce  these  conditions  all  that  is  re¬ 
quired  is  to  allow  the  nodal  values  to  "float". 

For  the  engine  inlet  problem,  the  cross  section  of  the  model  is 
stibdlvlded  into  a  finite  number  of  triangular  elements.  R  typical 
section  is  illustrated  in  Figure  5.  The  assumption  is  made  that  the 
scalar  eizlmuthal  potentials  V2u:y  line^u:ly  over  each  of  the  elements. 

This  necessitates  the  use  of  linear  basis  functions  defined  on  each 
node  of  every  element. 


Figure  5.  Typical  Finite  Element  Section  of  Model  Cross  Section 

These  pyramidal  functions  will  have  a  value  of  unity  on  a  specified  node 
and  a  value  of  zero  on  neighboring  nodes.  For  the  i^^  node  on  the 
element,  the  linear  basis  function,  ,  is  defined  as 

♦lO'W  +  =l,t* 

The  real  constemts  a^  b^  and  c^  ^  for  i  ■■  1,2,3  on  each  element 
may  be  evaluated  by  means  of  the  following  equations: 


3 


(19) 


(20) 

(21) 

(22) 

(23) 

(24) 

(25) 

(26) 

(27) 

(28) 


and  is  twice  the  area  of  a  triangle  (7:547) .  These  computations  assume 
a  cloclcwise  ordering  of  the  node  numbers  on  each  element. 

The  azimuthal  field  components  within  the  inlet  model  may  now  be 
expressed  as 

-  _  N 

4  •  e^(R,Z)  -  a^4j^(R,2)  (29) 


14 


♦  .  h^(R,Z)  -  a^*^(R,2). 


ccnplex  nodal  values  of  ^  and  ^  at  a  particular  node,  1,  are 
represented  by  euid  respectively  (6:204). 

Ibe  finite  element  approach  begins  with  the  substitution  of  the 
asimuthal  field  expeuislons  In  Bqs  (29)  emd  (30)  Into  the  Integral 
functional 

P  -  f  (31) 

8 

idtere  L(R,Z,4)j^,i|i2/V^^,V^2)  Is  defined  by  Eq  (14).  This  results  in 

N  N 

'  -  /'n  “r  ''•“Jl’fiV  ♦ 


«  N  N 

n  ♦  X  7  ♦  ^^^iSl^i^i^  * 


n  II 

^^r  ^^^lSl®i*i^  -  n  J  X  7(RJ^Yi*i) 


®r  *^^i-l''^i*i^^  -  ^i-l®i*i^^ 

The  stationary  point  of  F  Is  obtained  by  differentiating  the  functional 
with  respect  to  each  of  the  nodal  coefficients  and  setting  it  equal  to 
sero. 
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Therefore, 
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for  n  ■  1/2,  .../  N.  For  the  inlet,  only  the  case  for  n=l  will  be  examined 
In  addition,  it  will  be  assvmed  that  p^.  «  e^.  •  !•  Equations  (33)  and  (34) 
then  became 


r*i  R  •  V(R4J  -  R*,*_  dRdZ 
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(36) 


Equations  (35)  emd  (36)  cure  Integrated  only  over  the  elements  connected 

to  node  m  emd  relate  v  ^uld  3  to  adjacent  nodal  values  of  y.  and  6. . 

mm  11 

For  the  mesh  Illustrated  in  Figure  5,  14  adjacent  nodal  values  will  be 
related  to  one  ^mother.  These  equations  will  produce  a  2N  x  2N  sparse, 
banded  symmetric  matrix,  which  is  characteristic  of  cuiy  finite  element 
computations . 

All  integrals  to  be  evaluated  in  Eqs  (35)  and  (36)  are  of  the  gen¬ 
eral  form 


P^g  -  /  dRdZ 
Q  .  f  ^  dRdZ 

«r8  J  RiS-l 


(37) 


(38) 


A  simple  method  for  computing  both  types  of  these  integrals  exists  by 
invoking  the  two-dimensional  Stoke' s  theorem  «uid  converting  the  surface 
Integrals  into  line  integrals.  The  general  formula  is  given  by 

J  g(R)Z®  dRdZ  «  f  g(R)Z^*^lR)  dR  (39) 

8  C 

where  the  integration  is  performed  in  a  clockwise  maumer  around  the  tri¬ 
angular  elements.  The  factor  g(R)  is  any  £u:bitrary  function  of  R  (6:205). 
If  Z  is  expanded  in  terms  of  R  l.e.  Z  >■  mR  b  where  m  ■  slope  and  b  »  Z 
intercept,  the  form  of  Q^g  is 
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(R  -1) (R.+i)  ^  n 

-StoC  ?  .,,A")  -  2  ^ 

(R^-1) (R^+1) 


Higher  order  integrals  of  this  form  may  be  generated  by  the  recursion 
formula 


a 


R  o  R  -  R 


(6:205). 

In  general,  the  homogenous  system  of  linear  equations  defined  by 
(35)  euid  (36)  does  not  have  a  unique  solution.  However,  some  values  of 
nodal  coefficients  along  the  exterior  boundcuries  of  the  inlet  model  will 
be  specified  as  Dirlchlet  boundary  conditions.  These  boundary  conditions 
can  be  used  to  render  a  unique  solution  to  the  syst^  of  equations.  Before 
the  enforcement  of  boundary  conditions,  the  form  of  the  linear  system  is 


■^,2 

■'1,3 

•  ^,2N 

^2,1 

2,2 

■'2,3 

•••  *2.2H 

*^3,1 

■'3.2 

■S.3 

•••  *'3,2H 

^,1  *^2N,2  *^2N,3 


*^2N,2n|| 


\Aiere  the  K's  are  generated  by  integration  of  the  Idgranglem  over  the 
elements  £uid  the  y's  and  3's  are  the  values  of  the  nodal  coefficients. 
Mow  if,  for  example,  3^  is  specified  as  a  Dirlchlet  boundary  condition, 
and  has  a  value  of  V^,  the  system  may  be  modified  to 
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Figure  6.  Solution  Domain  for  Engine  Inlet  Model 

Boundary  represents  the  opening  of  the  inlet  model.  Along 
this  segment,  Dirichlet  conditions  for  the  fields  incident  on  the  mouth 
of  the  duct  cure  specified.  Due  to  the  cylindrical  nature  of  this  prob¬ 
lem,  these  conditions  are  in  terms  of  the  Bessel  functions  J*  (K  R)  for 

n  c 

e.  and  J^{K  R)  for  h.. 

The  outer  wall  of  the  waveguide  is  labelled  (§)  .  Because  of  the 
perfectly  conducting  nature  of  this  surface,  the  t^ulgential  coiq>onents 
of  the  electric  field,  e^  and  e^,  must  be  identically  zero.  In  terms 
of  aalmuthal  field  quantities  only,  the  equivalent  bo\md£u:y  conditions 


e^  -  0  (46) 

and  •  n  ■  *»  0.  (47) 

The  bladed  structure  of  the  ferns  in  the  initial  stages  of  a  com¬ 
pressor,  is  approximated  by  the  back  plate,  ©.  As  stated  eeurller, 

% 

azimuthal  current  flow  on  the  face  of  this  plate  is  to  be  eliminated. 
This  requires  enforcement  of  the  boundeury  condition 
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cos  6  -  sin  6  »  0 


(48) 


or  equivalently 

(■^  -  R  ■^)  cos  S  -  (1^  +  R  i||t)  sin  6.  (49) 

Additionally,  it  is  required  that 


e_  cos  £  =  e  sin  5  ■>  0 
R  Z 

vdiich  implies 


sin  £ 


_  SRe^  _  _ 
SR  ^ 


cos  £ 


(50) 


(51) 


Since  the  cone,  (§) ,  like  the  wall  of  the  waveguide, ,  is 
perfectly  conducting,  it  hats  similar  boundary  conditions.  Ihese  con- 
ditlcsis  may  be  e}q>ressed  as: 

e^  -  0  (52) 

and 

cos  »  (  -  -  R  ^^)  sin  (53) 


for  0  <  a  <  ir.  Ata*0  and  a  «  ir  the  boundaury  conditions  axe 

identical  to  those  for  (f)  which  aure  given  by  Eqs  (15) ,  (16) ,  auid  (17) . 

Onimoment  Method  Formulation.  The  general  process  of  applying  the 
unimoment  method,  consists  of  several  steps.  First,  a  form  of  the  field 
variations  is  assumed  along  each  exterior  boundary  of  the  inlet  model. 
Then,  using  these  Dlrlchlet  boundary  conditions,  the  resulting  Interior 
fields  aire  computed  using  the  finite  element  implementation  of  the  CAP 
equations.  This  process  is  repeated  N  times  for  N  different  forms  of  the 
field  variations.  From  these  N  solutions  a  set  of  weighting  functions 
are  calatulated  to  enforce  the  Neumauui  boundary  conditions  in  the  mean 


sense 
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AS  an  example,  consider  the  boundary  conditions  along  the  exterior  wall 
of  the  circular  cylindrical  waveguide  given  by 


Assume  that  h*  varies  as  sin  along  this  surface,  for  n  ^  1,2,  ...,M 
^  Zo 

where  N  is  the  total  number  of  nodes  along  the  boundaries  of  the  model. 
Appropriate  field  variations  along  the  other  boundairies  must  be  also 
enforced  simultaneously  with  this  one.  For  each  n  a  solution  is  obtained, 
thus  specifying  0^,  a^,  b^  and  c^  for  each  node.  For  the  total  field 
solution,  which  is  a  linear  combination  of  each  of  the  n  solutions,  it 


is  required  that 


I  If  .  0 

nSl  n  3R 


along  the  wall  of  the  waveguide.  Ihe  W^'s  are  the  weighting  functions 
for  the  n^  solution.  These  equations  are  than  enforced  using  testing 
functions  such  that 


(^1  az-0 


and  consequently  enforce  the  Neimann  boundary  condition  in  the  mean 
sense.  These  equations  will  generate  a  N  x  N  matrix,  and  produce  a 
linear  system  of  equations 


B  W  -  0 
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To  produce  a  unique  solution,  the  right  heuid  side  of  the  equation  must 
be  non-zero.  For  this  reason  some  small  6  must  be  used,  6  <  <  1,  so  that 

B  W  «  £  .  (59) 

This  is  the  general  method  to  be  used  to  enforce  Neumcuin  boundary  condi¬ 
tions. 
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Ill .  Results 


The  results  presented  in  this  section  concern  only  the  validation 
of  the  software  developed  for  the  finite  element  implementation  of  the 
coupled  azimuthal  potential  differential  equations.  Appendix  A  contains 
a  listing  of  the  program  used  for  this  procedure.  In  addition,  a 
description  of  all  variables  used  in  the  program  is  given  in  this  appendix 
To  Insure  the  program  was  correctly  ccmputing  the  azimuthal  field 
cco^nents,  it  was  necessary  to  enforce  boundaury  conditions  for  a  known 
wave  solution.  For  this  reason,  boundaury  conditions  for  a  standing  wave 
in  a  circulau:  cylindrical  waveguide -terminated  by  a  flat  conducting  plate 
were  used.  This  required  a  minor  modification  of  the  inlet  model.  The 
angles  o  and  5,  defined  in  Figure  6,  were  set  to  zero  emd  conducting 
boundaury  conditions  estadillshed  along  the  back  plate.  As  stated  eaurlier 
the  boundary  conditions  for  a  conductor  are 

e^  ■  0  (60) 

and 

V(Rh^)  •  n  -  0  (61) 

Since  it  is  not  possible  to  directly  enforce  the  Neumaum  condition  given 
Equation  (61)  with  the  minimized  Lagranglan  functional,  the  Dlrchlet 
boundaury  conditions  for  a  standing  wave  were  applied  along  all  conducting 
surfaces . 

For  a  standing  wave  in  a  clrculaur  cylindrical  waveguide,  the  elec¬ 
tric  auid  magnetic  fields  must  satisfy  the  following  set  of  equations: 

-  A  P/a)  sin  Kg  Z  e^*  (62) 


(63) 


64 


1  a^'Hz 

1  1  3^Hz 

^  a(^3Z  (9:415) 


Only  azimuthal  field  quemtlties  are  being  calculated,  so  it  is  rally  neces¬ 
sary  to  determine  the  forms  of  E.  and  H..  Using  the  2dx>v’e  relations  it 

9  9 

is  easily  shown  that 


^  jwy  p/a)  sin  KgZ  e^* 


For  all  computations,  the  azimuthal  dependence  is  suppressed  and  A  is 
assumed  to  be  one. 

As  a  first  test  case,  a  section  of  waveguide  was  evaluated  with  a 
radius  of  .SXg  and  a  length  of  l.SXg.  The  finite  element  division  con¬ 
sisted  of  96  trleuigles  with  65  nodes.  This  resulted  in  a  trlangulau: 

2 

■esh  density  of  126  triangles  per  .  Plots  of  the  on-axls  results 

for  the  magnitude  euid  phase  of  e.  and  h.  are  presented  in  Figures  7,  8, 

9  9 

9,  and  10.  In  general,  as  is  quite  evident,  the  results  ue  not  promis¬ 
ing.  As  a  measure  of  the  error,  the  squ2u:e  root  of  the  average  squared 
error,  i.e. VsCx-x)^/^,  was  evaluated  for  the  magnitudes  of  both 

azimuthal  field  quantities.  This  resulted  in  a  value  of  .462  for  e.  and 

9 

.524  for  h .  (all  values  were  normalized  to  one) . 

9 

It  is  believed  the  results  frran  this  first  trial  are  so  poor  because 

the  trlangul^u:  mesh  was  not  sufficiently  dense.  Morg^ul  states  that  a 

2 

mesh  density  on  the  order  of  800  triangles  per  X^  will  produce  an  error  of 
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Magnitude  of  (along  axis) 


Figure  7.  Comparison  of  Theoretical  and  Computed  Results 
for  Magnitude  of  e^i  Low  Density  Mesh 


Phase  of  (along  axis) 


about  2%  (6:209).  Consequently,  as  a  second  trial  case,  the  mesh 
density  was  increased  to  this  value.  The  radius  of  the  waveguide  section 
evaliiated  was  .31^  with  a  length  of  .41^.  On-axis  results  are  presented 
in  Figures  11,  12,  13  cuid  14. 

In  the  sense  of  the  error  measurement  being  used,  these  results 
2u:e  significantly  better  than  those  for  the  low  mesh  density.  The  sc[uare 
root  of  the  average  squared  error  produced  values  of  .137  for  the  magni¬ 
tude  of  e.  and  .223  for  h,.  In  addition,  the  phases  of  both  e.  and  h. 

♦  ♦  <p  ♦ 

eu:e  very  close  to  the  theoretical  values. 

As  a  final  trial,  the  length  of  the  waveguide  section  was  decreased 

to  .2Xq  with  the  same  radius  as  the  leist  case.  With  96  triangular  ele- 

2 

ments,  this  configuration  had  a  mesh  density  of  1600/Xo  *  Figures  15, 

16,  17  and  18  contain  the  results. 

Again,  in  the  sense  of  the  error  measurement,  the  results  continued 

to  improve.  The  average  squared  error  for  e.  is  .0334  zmd  .108  for  h.. 

♦  ♦ 

For  both  azimuthal  fields  the  phases  are  very  good  approximations  to  the 
theoretical  results. 

Since  none  of  the  trials  generated  a  good  approximation  to  the  magni¬ 
tude  of  the  standing  wave  in  the  interior  of  the  model,  the  wave  equation 
for  an  arbitrary  node  was  evaluated.  The  Interior  of  the  waveguide  is 
source  free  in  which  case  the  wave  equation  is 

V^E  +  =  0  (69) 

Expanding  in  terms  of  e .  and  h .  this  bec<»ies 
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Figure  13.  Comparison  of  Theoretical  and  Computed  Results 
for  Magnitude  of  h(h;  Medium  Density  Mesh 


phase  of  (along  eucls) 


Figure  14.  Compwison  of  Theoretical  and  Computed  Results 
for  Phase  of  h^i  Medium  Density  Mesh 


Magnitude  of  (along  axis) 


Figure  15.  Conparison  of  Theoretical  and  Computed  Results 
for  Magnitude  of  e^i  High  Density  Mesh 


Figure  18.  Comparison  of  Theoretical  and  Computed  Results 
*  for  Phase  of  h^;  High  Density  Mesh  * 


Since  a  linear  approximation  to  the  fields  Is  being  used,  the  seccxid 
derivatives  do  not  exist.  For  this  reason,  the  assumption  was  made  that 
the  fields  were  smoothly  varying  and  central  differences  were  used  to 
calculate  e^,  the  results  was  .1059.  This  indicates  the  wave  equation 
is  approximately  satisfied  and  the  program  is  functioning  properly. 


0 
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IV.  Conclusions  cuid  Recommendations 


Conclusions 

The  algorithm  used  for  the  finite  element  implementation  of  the 
Cap  differential  equations  functioned  pr(^rly.  For  the  three  mesh 
densities  used,  the  error  measurement  indicates  the  computed  results 
2u:e  converging  to  the  theoretical  values  as  the  mesh  density  is  in¬ 
creased.  A  logarithmic  plot  of  these  errors  is  given  in  Figure  19. 


Error 


Mesh  Density 


Figure  19.  Plot  of  the  Error  Measurement  for 
Increasing  Mesh  Densities 

This  plot  indicates  that  an  extraordinarily  dense  mesh  {  7500)  would 

be  reuqlred  to  achieve  acceptable  results  (2%  error) .  Fpr  an  inlet 
model  with  reasonable  dimensions ,  this  mesh  requirement  is  prohibitively 
luge  from  a  computational  standpoint. 


the  most  likely  cause  of  this  problem  is  the  enforcement  of  Dirichlet 
boundary  conditions  for  the  magnetic  field  along  a  conducting  surface 
Instead  of  the  Newnuum  condition.  However,  the  Lagrangian  functional  for 
the  CAP  equations  caumot  accommodate  the  Neumann  boundary  conditions. 
Therefore  the  finite  element  calculation  of  the  coupled  azimuthal  potentials 
in  the  engine  inlet  model  must  be  ckbandoned. 

BecCTanendations 

The  finite  difference  method  could  be  used  to  Implement  the  CAP  dif¬ 
ferential  ecpiations  directly.  Additionally,  the  Neumeum  boundary  condi¬ 
tions  for  the  magnetic  field  can  be  enforced  quite  readily.  This 
numerical  method  should  be  Investigated  extensively  as  a  replacement  for 
the  finite  element  method. 
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APPENDIX  A 


This  appendix  contains  a  listing  of  the  computer  program  used  for 
in^lementing  the  finite  element  method.  Definitions  of  all  major  var¬ 
iables  used  in  the  program  £u:e  provided. 


Veuriables 

ROORD(l)  -  Real;  contains  the  R  coordinate  of  the  Ith  node. 

ZCORD(I)  -  Real;  cont2d.ns  the  Z  coordinate  of  the  Ith  node. 

PLYNMIi  (I,J,K)  -  Real;  contains  the  K  coefficients  of  the  linear 
polynomial  of  node  J  c»i  triangle  K. 

Gammas (I)  -  complex;  contains  the  values  of  the  boundary  conditions 
for  the  azimuthal  electric  field. 

BETAS (I)  -  Complex;  contains  the  values  of  the  boundary  conditions 
for  the  axlmuthal  magnetic  field. 

GAMNUM(I)  -  Real;  contains  the  node  number  of  the  boundary  condi- 
tlcms  specified  in  GAIQIAS. 

BBTNUM(I)  -  Real;  contains  the  node  nuhmer  of  the  boundary  condi¬ 
tion  specified  in  BETAS. 

EQOTN1(I,J)  -  Complex;  contains  all  equations  minimized  with 
respect  to  the  electric  field  nodal  values. 

BQUTN2(I,J)  -  Complex;  contains  all  equations  minimized  with 
respect  to  the  magnetic  field  nodal  values. 

ANSNR(I)  -  Complex,  contains  the  final  computed  values  of  the 
azimuthal  fields. 

NUMEli-lnteger;  the  number  of  elements  in  the  system. 

IfOMNOD- Integer;  the  number  of  nodes  in  the  system. 

M(X}BS(I,J)  '  Integer;  contains  the  node  numbers  on  each  triangle. 

INTGRLd)  -  Real;  contains  the  results  of  all  integrals  of  the 
form  R^Z^. 

RSX^(I)  -  Complex;  contains  the  results  of  all  Integrals  of  the 

form  Sjli  . 

R2-1 

SYSTEM (I, J)  -  Complex;  contains  the  system  of  equations  defined 
in  the  Lagrangian,  before  inversion. 
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